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ABSTRACT 


Wing rock at high angle of attack is an oscillatory 
lateral-directional motion phenomenon known to exist in 
Sememot today's high performance tactical aircraft. The 
motion has been consistently characterized as a lightly 
damped Dutch-Roll oscillation attributable to asymmetric 
wing stall. However, evidence gathered from wind tunnel 
Simulations and at least one British study indicate that 
aerodynamic nonlinearities may be the source of wing rock. 
Regardless of the actual cause of the phenomenon, a study 
of wing rock has positive ramifications with respect to 
gaining a clearer understanding of the aerodynamics associ- 
ated with high angle of attack flight. This report presents 
the results of an investigation of wing rock which centered 
on the premise that two distinct nonlinear aerodynamic mech- 
anisms (aerodynamic hysteresis and a cubic nonlinearity in 
yawing moment) not only can cause wing rock but may drive 


Mamvora tinigecvcle ostillation as well. 
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LiSr OF SYMBOLS 


Wing span 


Dimensionless rolling moment coefficient (Ret) / dob 
Dimensionless yawing moment coefficient = (Y.M.)/qSb 


Mass moment of inertia. Subscript x, y or z denotes 
reference body axis. 


(Rolling Moment) /I,, pOsitive tends to lower right 
wing. 


es Moment)/I., positive tends to yaw right wing 
elit Ee 


Pollvgase = @P/dt (or? ) 
Dynamic pressure, = Pu’. 
Yaw rate = dP/dt (ory) 

Wing area 

Period of motion(2%/w) 


Freestream velocity 


Greek Symbols: 


fi 


| 


ai, Ss 6 
tl 


Sideslip angle, positive for wind vector coming from 
the right. 


Bank angle, positive for right wing down. 

Yaw angle, positive for nose moving to the right. 
Phase angle 

Dimensionless damping ratio 

Dimensionless nonlinearity variable 

ico peanplitude term 

Pigenwalwe sr Gatio of /- as noted in text usage 


Atmospheric density 





da = Aileron deflection 


A()#= Incremental change of a variable 
pee— Carcular frequency 
(w= undamped natural frequency) 
(wi = Yl- Ge a = damped natural frequency) 
Bebpscripts: 
( Jy = d( )48 .; partial derivative with respect to sideslip 


ane Le 


C97 90/8 Pdr oy 


Time Derivative: 


m= dd j/dt 
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I. INTRODUCTION 


Wamemrock at moderate to high angle of attack is a char- 
acteristic known to exist in some of today's high performance 
Mm@eerart, particularly those possessing the characteristics 
Seeethin swept wings and a slender fuselage. The terminology 
“wing rock" has been used to characterize an oscillatory 
Pateral-directional motion exhibited by some aircraft when 
operating near stall, or when flow separation is evident by 
airframe buffeting. Although the motion has been observed 
and documented extensively during flight tests, its exact 
cause has remained essentially unknown. In an attempt to 
explain wing rock, fleet aviators, test pilots and aerodynam- 
icists have consistently advanced the theory which has gener- 
ally concluded that the apparent random rolling moments 
manifested by alternating wing rock motion at or near stall 
are due to asymmetric spanwise loading, i.e., alternating 
Small Gondition along the wing span. It is not the intent 
of this thesis to dispute such arguments but rather to present 
new arguments based on the premise that aerodynamic nonlinear- 
ities may not only cause wing rock but may drive it to a limit 


@ycle oscillation. 


pe enol EXPERTENCE WITH WING ROCK 
Micmalth@@elas Nad experience as a test pilot at the Naval 
PetelesrenGenter and has personally observed what could be 


essentially described as a wing rock limit cycle while 
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Gomewetimme high angie of attack flight tests with the Northrop 
T-38A airplane. The T-58A is control limited longitudinally 
to the extent that it can not be made to achieve the so-called 
iersscic aerodynamic stall", i.e., complete flow separation 
mone theswing span. As a mesult, the normal l-g stall is 
defined when a maximum angle of attack is reached due to 
control-limiting, moderate-to-heavy airframe buffet and a wing 
rock motion are evident with an average amplitude excursion in 
bemk angie of + 20 degrees. This lateral motion occurs with 
essentially no heading deviation and will continue indefinite- 
ly, albeit accompanied by a large sink rate, until the angle 
Seeeateack is reduced. The buffeting and wing rock will 
immediately cease and a complete recovery can be effected 

with very little additional loss of altitude, when afterburner 
Ghrust is applied. 

The question of whether wing rock motion in the T-38 (F-5 
series probably behaves the same) is indeed a limit cycle 
oscillation caused by aerodynamic phenomena other than asym- 
metric spanwise loading has never really been seriously con- 
Sidered from a flight test standpoint. Instead, most test 
pilots and engineers associated with experimental/developmental 
flight testing quickly assume that an asymmetric stall condi- 
Promeis responsible, since the motion is typically accompanied 
Meech aGdeiOneas eyidenced by airframe buffeting. Most swept 
wing aircraft exhibit wing rock motion to one degree or another 
near stall, although few can actually be said to exhibit motion 


Dee eeomedomae iinet Cycle, Since mOSt aircraft are not control 


ia 


Mmime@ecmmtne dpproach to stall will often continue past that 
Pommeeneme an Oscillation similar to that exhibited by the 
iemommoht occur. The Vought A-/ and F-8 series are notable 
examples of aircraft that can quickly reach a point where 
@arectional stability is lost, resulting in an almost immedi- 
@ee departure from controlled flight, if the pilot does not 
respond quickly enough to reduce the angle of attack. 

Informal liaison with NASA Langley revealed that free 
imemont: model testing of the I[-38A had on at least one occasion 
confirmed the existence of a wing rock limit cycle. Although 
credibility of this report would have been enhanced had it 
been possible to conduct an analysis of lateral directional 
MetiOn at high angle of attack on a current production air- 
@ratt such as the 7-38 or F5, the necessary information needed 
to conduct such an analysis was not forthcoming owing to 
Peeprietary considerations. Thus, it was necessary to inves- 
Meeatc and illustrate the concept of a wing rock limit cycle 
using the textbook airplane, F-94, as presented in Blakelock 
(Ref. 1]. During the course of the research the stability 
derivatives of the production F-94 were modified as necessary 
to simulate desired test conditions. To the extent that the 
Povemcemiccomds a vehicle with which to demonstrate the concept 
of wing rock as a limit cycle, it was a satisfactory substi- 


Piecminm Lieu or a current production aircraft. 


B. A NEW THEORY PROPOSED 
Explaining wing rock as an aerodynamic phenomenon apart 


from the theory of asymmetric wing stall poses some difficulty 


De 
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Pionioeerom tne Lact that the motion closely resembles a 
lightly damped Dutch-Roll oscillation. Evidence obtained 
meemeclosely controlled wind tunnel experiments [Ref. 2] and 
Seeervattons made during flight tests seem to suggest the 
Deesibility of a limit cycle oscillation wherein the amplitude 
and period of the motion may result solely from nonlinear 
aerodynamics. This contrasts with the response of a lightly 
damped Dutch-Roll oscillation where the amplitudes are a 
mametion Of the initial conditions. 

Some difficulty arises when attempting to identify an 
aerodynamic nonlinearity that may contribute to wing rock 
momron in wind-tunnel or flight tests since most dynamic 
Bemirysis techniques, including aircraft parameter identifica- 
mem, are Oriented towards a linearized model of the aircraft. 
Presumably, a nonlinear model could be assumed for parameter 
identification purposes but a question of uniqueness would 
arise since several completely different plant models might 
produce equally valid time history matches with actual flight 
feweae Therefore, a significant step at this time would be 
mremrceoemwition OL wing rock motion as a limit cycle oscilla- 
Gion induced by aerodynamic nonlinearities. 

POssei net. 5], has Shown that a wing rock limit cycle can 
SseemmicieanOmlinearity of cubic form is present in either 
Vopiowon rolling moment with respect to sideslip angle. 
BOaomuscamtme Kryloft-Bogoliubot method (KR method) with the 
assumption of slowly varying parameters to analyze the mathe- 


matical model of the motion described by the test airplane. 


lS 





Pepommuen OL tChewanalysis presented herein explores further 
the work done by Ross with the objective of obtaining a 
Samplified form of the theory she proposed as well as ampli- 
meee Upom it. Preceding this, however, will be an analysis 
and discussion designed to support the claim that a wing rock 
emat cycle may also be the result of a nonlinear mechanism 


referred to as aerodynamic hysteresis. 


ee MOTIVATION 

MnewsStuady Of Wing rock and its causes has direct applica- 
tion to an analysis of stall and post-stall gyrations which 
have long been difficult to predict and control on many 
tactical aircraft. The Vought F-8 and A-7 are notable exam- 
ples. Increased knowledge in this area would surely enhance 
the ability to predict aerodynamic behavior, thus enabling .- 
an earlier detection of undesirable flying qualities, or 
design improvement or, at least, a minimization of mission 


limiting or annoying characteristics. 
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II. METHOD OF INVESTIGATION AND EQUIPMENT USED 


iime Imvestigation Of wing rock involved the use of the 
Peewee (Electronics Associates Incorporated) model 580 analog 
computer as well as the Naval Postgraduate School's Depart- 
femme of Aeronautics digital computing facilities, namely, 
the Hewlett Packard, HP-9830. Mathematical models represent- 
ing the pertinent equations of motion were programmed on, and 
meaedlts Compared from, both machines to ensure validity of 
Paocramming. Ihe output from these computing machines con- 
sisted of time histories of the motion exhibited by the test 
aircraft (F-94) when released from some initial condition. 
Semele data were recordéd by a graphical plotter or strip 
@m@art recorder attached as a peripheral device to the main 
Mamiputing Units. These time history plots were used to make 
been GUualitative and quantitative assessements of the motion, 
including measurements of critical damping ratio, frequency 
and/or period. Digital as well as graphical output was ob- 
tained using the HP-9830. Time history plots appearing herein 
were initially made using the analog computer, but final copies 
Pf oeemoemammeds strom the HP-9550 plotter as it produced a more 
compatible form for this report. The analog computer was 
Sumgeortor to the digital machine in terms of its ease of pro- 
gramming or modification of the program during execution where 
Pe imp eocmpucecO vVicw readily the effects of altering certain 


LVetirohammparbameters . 
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ace 
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The mathematical model used for laboratory tests consisted 
of the linearized lateral-directional equations of motion with 
certain simplifications to be presented later. As mentioned 
earlier, the F-94 airplane as presented in detail in Blakelock 
meen. 1) was used as the test vehicle for conducting laboratory 
eemulations of wing rock motion. In its production configura- 
jen the airplane exhibits no inherent instabilities with 
respect to the lateral-directional modes of motion. Modifica- 
tions were made to its production stability derivatives as 


necessary to conduct desired simulations. 


16 


1948 





[11]. RESULTS AND DISCUSSION 


The discussion will focus on the particular lateral- 
directional motion exhibited by some aircraft within their 
meepective flight envelopes and referred to hereafter as wing 
meek. The first part of this section will be devoted to a 
presentation of the analysis and investigation of wing rock 
Semcaetimit cycle induced by aerodynamic hysteresis. The 
second part will present a similar analysis and investigation 
cme rock as a Jaimit cycle induced by a natural instability 
immene aircratt that is constrained by a nonlinearity in one 
Or more stability derivatives. The reader is advised to refer 
to pages 8 and 9 for definitions of symbology used throughout 


mae 2Ollowing discussion. 


feet GENERAL LIMIT CYCLE OSCILLATION 
wom weve lesosciliation is ssimilarn physically to the 
feemen described by the typical undamped harmonic oscillator 
mepuesented by: 
ne ae 
Wien has solution: 
x(t) = As (vx t +) 
x(t) = A Ye Cos (Yr t + B) 


The motion amplitude of the linear undamped system will depend 
UpOMmumemiInatial Conditions, and hence different initial con- 


ditions will result in different amplitude histories. 


ile 
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=a = ite. 
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The linear undamped system will, of course, exhibit no damping 
Delavaor or the motion. By contrast, a nonlinear system that 
dememsenaces a limit cycle behavior will have the same ampli- 
fede description for motion regardless of the initial condi- 
tions. A system that is weakly nonlinear with a limit cycle 
behavior will appear to have the properties of a linear system, 
however, the characteristic of motion amplitude regardless of 
mercial Condition is a distinctly different attribute of the 
nonlinear system. One should recognize that the above differ- 
ences are somewhat heuristic and a mathematician would char- 


acterize the distinctions with considerably more detail. 


Pee ee ERODYNAMIC HYSTERESIS 

Wind-tunnel tests on aircraft models often show that the 
rolling moments measured on a symmetric model, when yaw angle 
memeoeno and angle of attack is increased, are approximately 
Bomueovalued until the angle of attack approaches close to 
cea. ROlling moments tend to be erratic at or near stall, 
and it is customary to attribute these traits to slight 
asymmetries in the model or tunnel flow, with a resulting 
early stall of one wing relative to the other. A recent wind- 
tunnel test [Ref. 2] upon a sting-mounted model having a strain 
gage balance with appropriate electronic filtering, used the 
test technique of conducting a slow, continuous yaw angle 
Vcc eeluelieteanconplecte cycle of yaw motion. When the cyclic 
yaw sweep was conducted at or near stall angle, both rolling 
and yawing moments showed incremental shifts in value that 


appeared to be related to the direction of the slow yaw sweep. 
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Tryeseiueeemon rOlling moment, when plotted versus yaw angle, 
gave evidence of aerodynamic hysteresis. 

Sumely Stated, aerodynamic hysteresis is a type of behavior 
Waerein an aircraft's response to a given input is such that 
it does not follow the same path returning to equilibrium as 
meedid when initially displaced. Of importance here is the 
immer that the aircraft remains positively stable. Figure l 
[emai tllustration Of stall hysteresis wherein the airplane 
peepetriences a sudden decrease in lift coefficient as a result 
Sescomplete flow separation along the wing span at stall angle 
of attack. The airplane recovers to a lower lift coefficient 
ummole OL attack Only to repeat the cycle if the controls 
ape held in a pro-stall condition. This type of motion is 
Gruen cererred to as “stall porpoising' and is common to most 
m@eeratrt With high aspect ratio wings. This example of longi- 
miemnial hysteresis 1S presented to the reader as an aid in 
Mueterstanding the discussion to follow pertaining to lateral- 


@rrectional hysteresis. 






| a 


a LONG/TUDINAL 
HYSTERESIS LOOP 


LIFT COEFFIC/ ENT ~C, 


ANGLE of AITACK~ & 


Figure 1. Stall porpoising and longitudinal 
hysteresis. 
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(ieeminaiysisspresented here will include some simpli- 
fications where deemed appropriate. The interest is to use 
State variable format to represent the mathematical model of 
the pertinent lateral-directional motion. The state variable 
format will be extended in an illustrative sense to demonstrate 
a wing rock limit cycle due to aerodynamic hysteresis. It is 
assumed that the lateral-directional motion of the aircraft 
occurs in such a manner that the airplane remains on a straight 


flight path, hence yaw and sideslip angles are related by: 


ee an 


Further, it is assumed that the aircraft 1s operating at a 
fixed angle of attack with all control surfaces fixed. 


Therefore, 


Ala= > yawing moment 


_b 
- g 5blC,, 8 + 2T(G2 Cre?) S (02a) 
and 
pt = = rolling moment 
= ¢ ShLGy, 8 Lr (G,.2 - GA) 1 (02b) 
Since #2 = ge = - gf. = a and R = ig equations (02a) 
and (02b) can be simplified to: 
B= -Wht*MS -WyeP (03a) 
and 
P= 4pfh KB -hoP (03b) 
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Further modifications have been made yielding a dimensional 
Pore tle lifcarizea lateral-directional equations of motion. 
Piewassumption that a straight flight path is maintained 
Merectively G€liminates any influence from sideforces and the 
usual accompanying spiral mode. Such an assumption is reason- 
able and consistent with flight test observations of wing rock 
motion. 

The freely vibrating system represented mathematically 


mm equations (03) may be cast into state variable format as: 


)x5 = [a] U4 


wren Corresponds to the shorthand notation for the following: 


a On ee fp 
6 = “Ng NV, Mo 3 (04) 
Ve La oh weal 


Pmearrad equation or identity 828 is necessary to accommo- 
date all the state variables in the matrix format. 
Romdmemess Druet ly, $x$ represents the state vector 
which contains the components of sideslip angle, #, sideslip 
velocity, %, and roll rate, 7. The plant matrix [aA] is of 
Order 3 and contains the dimensional stability derivatives 
(standard NASA convention) which resulted from the previous 


analysis. These derivatives are explicitly defined as follows: 


a, = Ar = O ; A,» = A 


Zr, = -Ng = 38 


\ 














b -¢I5d 
y= NK = ee Cy ea Cf 
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Equation (04) represents an eigenvalue problem of a freely 
Vibrating system in which the eigenvectors are the modes of 
motion and the eigenvalues contain the corresponding char- 
acteristics of damping, natural frequency, and period. In 

its present linear form the system described by equation (04) 
has a solution consisting of three roots (one real and one 
complex conjugate pair) corresponding to the roll convergence 
(subsidence) and Dutch-Roll modes respectively. Again, making 
use of the state variable format the solution to the Brae iur 


Vibrating system may be expressed as: 


fede [oe] fos} 


Pieme [oe (t)] is the state transition matrix and /x(0)f is 
the vector of values for the state variables at time equal to 
zero. In essence, the state transition matrix is used for 
modal analysis to project or propagate the aircraft's motion 
(WHecurecaity trom Some Starting values (initial condition) to 
another set of values at a specified time greater than zero. 

The equations of motion for the dynamic system with 
Control applied is given by: 

fx} = [a] {x} + [a jus (06) 


ee 





Wieme the Control matrix [8] contains the control power 
derivatives for yawing and rolling moments associated with 
the aileron and rudder controls. The vector ful contains 
Gieeme@enrtudes of the inputs for the respective controls. 


The general solution to equation (06) is expressed as: 


Fxcay - [ ore) {xc + [60] | [er] fo] fur} im (07) 


where the integral term 1s a convolution integral. This is a 
Memeral solution in the sense that it applies to any of the 
Seamaard anput methods, i.e., step, ramp and impulse. The 
amalysis and experimentation to which this section of theory 
meoletes Was restricted to the case of step inputs, which is 
probably the most common control input method used in aero- 
GymemicsS £Or Stldying aircraft dynamic stability. 
Uiiewspecit1e6 Solutron to equation (06) for a step 


input of control is expressed as: [Ref. 4] 


face} = [602] fueo} [a] [or 2] [a] fu} 8 


/ é . 
where [&] and [2%] are the plant inverse and identity matrices 
meopectively. 

Figure 2 contains a representative form of the aero- 


dynamic hysteresis assumed in the analysis of rolling moment. 


Zo 








Figure 2. Assumed form of roll hysteresis. 


The motion of the system corresponds to a periodic variation 
of the state variables such that the amplitude of the sideslip 


angle is assumed to be bounded by: 


£48428 


Referring to figure 2, the hysteresis loop is defined by the 
dashed lines. The value of the rolling moment corresponding 
to a particular sideslip angle will be the sum of the produc- 
mammaire ratrt’s linear roll-due-to sideslip and the additional 
hysteresis term, which has constant magnitude but varies in 
Seemmeaeecording to the sign of the sideslip velocity. This 
latter condition alludes to the manner in which the increment 
of rolling moment was applied, i.e., through a relay action 
Peoemmvion (+) driven by the sign of sideslip velocity. Side- 
Seapevclociwty awill naturally be plus or minus in sign according 
Omer dimectioOn Of yaw, hence the increment. of rolling moment 
was added or subtracted accordingly. Thus, the rolling moment 


due to sideslip as shown in Fig. 2 may be expressed as: 


Roll d@we= to sideslip = L (8) = Le f& + AL * SIGN (8) 


24 


— ea 





MeetucezZeschnoula not be construed to represent the only form 
meee roll hysteresis may take but was used in this analysis 
Gomeitetrate the Concept. Nothing has been stated explicitly 
faeinpliedsup tO this point with respect to the direction 
@emmmed by the hysteresis loop as the aircraft oscillates in 
lame i. e¢., Clockwise or counter-clockwise, but this will be 
meorrcicd Mater in the presentation of expemimental results. 
Although not shown, yaw moment hysteresis may be assumed to 
fomoe similar form as that presented in Fig. Z. 
[ImemctcwmeCmnnedi Ctethe behavior of the aircraft wt 
Wars MeCessary to solve equation (08) after some simplifications 
mem@iared to facilitate the ease of computation were implemented. 
Under the assumption of a limit cycle, the motion is periodic 


and the amplitude remains constant, which can be expressed 


fxcef = Jxce +7) 


meee le tsmcrne period, Referring back to the column vector 


mepart as: 


of state variables in equation (05) and with a convenient 


Smeace Of the Cime Origin: 
(aO- Gee a) -<wea= 2 


Rimemicmaenecessary and sufficient condition to assure proper 
establishment of the relay action on roll moment hysteresis. 


The forcing function portion of equation (06) may be expressed 


La] fuf a AL {3} 516 (8) (09) 


as. 


TAS 





where We = ae Ay 
Mime Value of unity in vector element u(3) of equation (09) is 
used to denote a step input of the roll moment increment, AL. 
This can be equated to a step input of aileron by the pilot, 
though roll hysteresis should not be considered as having that 
as its source, in keeping with the requirement for this analysis 
that controls be fixed. An external source May, Intact, be 
Mmeemorigin of roll hysteresis. Carrying the previous logic 
further, it is valid to consider the motion over a half-cycle 


Since periodicity once again requires that: 


Siw} 


7 X.(7/2) =0e 


fx(7%)} 


and X, (0) 


Ssupstituting these simplifications into equation (08) yields: 


fxcm)h = gerry f{xiot +[a] [acr)-z) [a] {4} AL 

ue] = 900m fea} =f" [sera] faa 
-[reocny fxio} =[4] [ecr)-z] [a] {af at 

{xcoy} = - [r+ sim] [aj [ar-2] [a] {ufAL (go) 


Peeteron (09) can be solved for the period T by numerically 
iterating for the value of T/2 which will cause X,(0), side- 
eee cuocity. CO Vanish. However, it is necessary first to 
solve for the state transition matrix [¢(+)] in order that 


all values are known except T/2. Coincidental to computing 
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the period, the maximum values of the state variables X4 (0) 
and Xz (0) (sideslip angle and roll rate, respectively) were 
determined, since these necessarily occurred at the state 
point where X, (0) woeoned. Computationally, this was a 
Sump ler method than integrating equation (07). A sample 
solution will be presented later. Upon examining equation 
eo, it should be noted that the period of the limit cycle 
in no way depends upon the sign of roll moment hysteresis; 
however, the amplitude of the state variables is a linear 
function of the magnitude of roll moment hysteresis, AL. 
The same procedure as outlined above could be applied to an 
analysis of yaw moment hysteresis. 
fee Analysis 

POtrereasons Piven eariiér in the introduction, the 
Peoreairyplane was chosen as the test vehicle for conducting 
Mperatory simulations of wing rock motion. In addition to 
being in the landing configuration, the following pertinent 


data relative to the F-94 was obtained from Blakelock [Ref. 


U = 135 mph = 198 £ps Ce, = .0278 RAD} 

q = 46.6 psf ay . 
"\ascaciaalilg (gy, = --0916 

S = 239 ft 

Peay. 5 ft Cy g = -105 

Ix = 7.16 x 10° slug-ft? Gy = -.055 " 
- 3 ge 

iher> 55.0] ae 10 a ie Con --219 " 

Cog = -,0487 RAD 

Gp =-.450 " 


For these values the plant matrix of equation (04) became: 


ea 


Ne: 


_— 





Q 1 0 


it 


[A] Zhi Sale 2491 0629 (10) 


mi Oe L ale 7 (tents Sha) i) 


a. Numerical Predictions 
The system of three first order, ordinary linear 
differential equations in (04) were solved by using the .com- 
puter program subroutine BASMAT, listed and described in Melsa 
and Jones [Ref. 5]. The Standard Fortran IV code was trans- 
lated into the BASIC program language code for adaptation of 
Peer) to the HP-9830 digital mini-computer. The eigenvalues 


of equation (04) using the plant defined by equation (10) were: 


DAs SECA 


> 
iI 

8) 
i 


hy 3 = b+ ic = -.12873 + 1.17554 SEC™* 


These values contain characteristic information pertaining to 
@aemewo modes Of motion which the aircraft exhibited as a 
freely vibrating system, i.e., no hysteresis imposed. These 
two modes were the roll convergence and Dutch-Roll modes. The 
roll convergence mode is indicated by the single real root, Ai> 
while the Dutch-Roll mode is represented by the remaining 
Sompilex roots, Ad 3: This information was needed for compari- 
son with the motion solutions obtained when hysteresis was 
applied (to be presented later). The Dutch-Roll mode indicated 


Biomol mis ecnaracteriStics with respect to period, frequency 


endsdamping of the stable motion. 
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KOs 0.12873 sec”! 


UW), = |}/- 4° Wry = 1.755 sae: 
Hence: 3S eles 
_ -l 
o8) = 1.183 sec 
n 
re goa 5 Sec. 


These calculations show that the Dutch-Roll mode is indeed 
stable, as expected for the production aircraft, with a damp- 
ing level and undamped natural frequency of 10.9% of critical 
and 1.1835 rad/sec, respectively. The roll convergence mode 
was not considered further, since the primary interest was in 
tine oscillatory mode Of MoOtren, 

Wiemstate transi tron Matrix, [@(7)], was obtained 


likewise from BASMAT and is presented below. 


De 
fece/ :[R]e** + Ie] COos(C%) + [6] sncee)/ =, (11) 
where 0.0114 ORO OZ G00 95 
[R] = ar 279 = 0 eau) 3.0) =O). Os 


10 Or SS Os) S) Lis 


0.9886 -0.0012 -0.0093 
[F] = 0.0279 1.0030 0.0228 


ale CL eo = eae: 0.0084 


0.1320 0.8531 0.0184 
Glew “ike VC eT) LOE 0.0086 


=O OF a0 Sale 5 O'S =a EZ or0) 





and aap, ¢ sderived trom Ay = a 


ms =—be + iC 


Peeme this Information input to equation (09), and applying an 
iterative solution technique, the period and amplitudes of the 
state variables for the cases of imposing roll and yaw moment 
hysteresis either independently or concurrently was obtained. 
The results for varying values of roll and yaw moment hysteresis 


@eempresented in Table I. 


TABLE I 


Resultesm@ot Numerical S@lution wt Eq. (09) for 
Different Magnitudes of Roll §& Yaw Hysteresis 


Type —_- Max Max Period 7° 


Yaw -4210 es LOZ eeeknon 
Rol] = gli) Shag See Oo 4 Do Ovan 
1 


Roll (om 
.3007 .0419 5.4219 


eet cc 03 
rer Oy Ie: 4 oboe OS! e554 12 1 


(1) Only condition for which roll rate,/, is not a 
maximum. All other values are necessarily maximums. 
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(te waterarive solution technique of interval 
ieeemioneas described by Gerald (Ref. 6] was used to solve 
meee (09) for the ™values in Table I. Two function values of 
Saeesite sign were initially obtained by inputing estimates 
@meene period. Thus started, the interval halving routine 
@emvetced in an iterative fashion to a solution. The computer 
peeeream used for this purpose and implemented on the HP-9830 
is presented in Appendix A. 

The state transition matrix, with its linear com- 
bination of the eigenvectors (mode shapes), allowed identifi- 
cation of the individual parameters comprising the Dutch-Roll 
mode. This identification follows from equation (11) and the 


following trigonometric identity: 


K(Bsnve+2 cos 6) = Kk sin (6+ 9) 


if) 


HAH SINO +4 Cos 4 
where 


K = A’ + B” = magnitude of motion vector in polar 


EOOndInates 


PD = Arctan B/A = phase angle of motion vector in 


DOmumCOOhCimiaces . 
Applying this mathematical form to the second term on the 
right of equation (11) results in: 


Lr], 


: Comet t) + Ie}: Si (C%) = K (2 SIN @ +h cos 8) K sin (a+) 


eZ 
where Mae = Gr + SIs) = magnitude of motion vector 
dp A aad h le of 
an . . = Arctan : = phase angle of motion 
od Chi vector 
ey = array locations. 


oul 





(eamowenis and appropriate initial conditions such that 


& (3) ie 
X (0), = ) ga) =) 
WP (0) 
the resulting relative magnitudes and phase angles of the 


elements comprising the total Dutch-Roll motion vector when 


normalized with respect to sideslip angle were: 


B& =1.00 fg= 0.0 Deg. 

| ; ff = 1.183 i ¢}. Pa = 96.25 Deg. 
2 Eee Py= 187.38 Deg. 
f =1.030 Pp = 91.13 


PeeiOuchmbanikeangle was not a state variable, its 
value can be determined as a function of time according to the 
@eeationship P(e) = [ve g/t , or approximated by B(t) = Ze 
The integration approach was used on both the analog and digital 
computers. Bank angle was introduced to gain a better under- 
standing of the motion's character, specifically with respect 
to the relative amplitude ratio of bank angle to sideslip angle. 
A positive value for the phase angle above indicates phase lead. 

b. Verification of Numerical Predictions 

The numerical computations presented in the fore- 
Perereeyiere Veriricd qualitatively by modeling the subject air- 
craft in the laboratory on the EAI-580 analog computer. The 
circuit required for the hysteresis model is shown in Appendix 
E. The same model using the system of equations given in 


equations (03) was programmed on the HP-9830 for the purpose 


ou 
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Seeomtatmineg Quantitative data. Digital propagation of air- 
craft motion was obtained by a numerical time-solution of the 
differential equations using the Euler Predictor-Corrector 
method given in Ref. [6]. The computer program used on the 
HP-9830 for this purpose is presented in Appendix B. A 
Separate program (Appendix C) for plotting output from the 
program in Appendix B was necessary to avoid overflowing 
available memory capacity due to large storage arrays used. 
Time history plots of aircraft response are pre- 
sented in figures 5 through 10. These computer generated 


plots were obtained by numerically integrating, as previously 


mentioned, the following governing equations of motion: 
P= hgh -L£r,G -4pP + Al & SIGN (8) 


These are equations (03) with roll and yaw hysteresis terms 
added. The third equation contained in matrix equation (04) 
was merely an identity and thus served no useful purpose in 
performing numerical integration. 

Tie weime history responses in tigures 3 through 6 
miGewanmranged to facilitate comparison of the characteristics 
of the two forms of motion, i.e., the damped Dutch-Roll and 
mepeemer cle OSCillation. The Dutch-Roll mode was not dupli- 
Gated) on all plots; therefore, it will be necessary for the 
reader to refer back to figures 3 through 6 for comparison 
purposes. All of the time history data were plotted with 
Sideslip angle as a reference parameter. 

Pigures 5 through 6 show the limit cycle resulting 


from the application of an increment of roll moment hysteresis 
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(AL=1) via the relay action mechanism. The relay action was 
Simulated on the analog computer by patching the output of 
sideslip velocity through a comparator and relay module whose 
Gmeput of + 10 volts was a function of the sign of #. ie 
same mechanism was simulated on the digital computer by making 
use of the predefined function, SGN, which returns a plus or 
minus 1 value according to the sign of its argument. 

Two rather unusual features that require some 
emempeOration are present in the time history plots. Specifi- 
cally, the unusual shape of the roll rate curve with roll 
hysteresis applied was believed due to the sudden application 
Of a rolling moment step input via the relay action. The 
lower of the two peaks per half cycle was judged the more 
accurate for analysis (refer figure 4). Secondly, a solution 
was obtained to equation (08) using a negative value of the 
mommies teresis, which except for a change in sign was identical 
to the solution obtained using a positive value of equal mag- 
Mimera@esrOll hysteresis (Table 1). However, both the analog 
and digital computers gave a time history response for the 
same value of hysteresis that was heavily damped or essentially 
ieamedisatter One cycle. The numerical prediction was there- 
tememassumed tO be in error, although this did not reflect on 
Siem bidity sol wtnewmathematreal modeling used. Rather, it 
MIomecnelided tiat the predicted solution failed to account 
-OMmeetenoveralieimtenaction of the production aircraft's sta- 
bility derivatives during a dynamic propagation of the motion. 


weed manor nemeway the hysteresis action would only traverse 
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imine COMmMeer-ClOCKWISe Girection about the loop in figure 3. 
A similar check was made on yaw hysteresis with identical 
ios tS. 

Figure 7 shows the limit cycle caused by applying 
yaw hysteresis via the same relay mechanism used for roll 
hysteresis. When compared to the corresponding time history 
plot for the Dutch-Roll oscillation in figures 4 and 5, the 
motion induced by yaw hysteresis in figure 7 bears a remarkable 
resemblance. The phase angle relationship and period (Table I) 
are identical to the Dutch-Roll. From this, one may conclude 
that yaw hysteresis is essentially a forced Dutch-Roll response. 

The combination of roll and yaw hysteresis reflects 
to some extent the relative influence of each on the total 
aircraft response. As one might intuitively expect, the result- 
ing motion was a compromise or average of the individual effects 
of roll and yaw hysteresis as previously seen. In Table I it 
was pointed out that the maximum value for roll rate was not 
necessarily achieved when sideslip velocity, 8, was zero. This 
was confirmed by the plot in figure 9. (Although not explicitly 
shown in figure 9, 6 was obviously zero when # was an extremum). 
There were many more combinations of roll and yaw hysteresis 
that could have been tried; however, this was not considered 
feasible nor productive. Again, the relay action can be seen 
in the behavior of the roll rate response in figure 9. The 
distortion of the roll rate response curve was due to the 
ScMotiecmetmects sot FOll and yaw hysteresis. Likewise, the 
phase shift for bank angle is a compromise between that seen 


in figures 7 and 5. 
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Roll and Yaw Hysteresis. 


42 


4H 


30 





Meth the exception of figures 6 and 8, the limit 
eceme Oscillations were not completely developed at the ter- 
fem@itroOn Of the time allotted for the numerical integration. 
irs was due primarily to the limited plotting space avail- 
able, whereas in figures 6 and 8 the phase plane plot was not 
so restricted, and clearly shows convergence or divergence to 
the limit cycle amplitudes after release from the initial 
conditions. 

For the sake of visualization, one could equate 
[mmeevalues Of roll or yaw hysteresis to a deflection of the 
MaesOpuiate COntrol surface, i.e., a nominal + 10% of avail- 
Mmencentrol power trom the aileron or rudder control surfaces. 
From the analysis conducted on the analog computer, with its 
real time ability to reflect changes in the problem variables, 
it became apparent after some experimentation that the cross 
coupling between roll rate and yawing moment through the 
adverse yaw stability derivative, VY, though weak, was respon- 
mmercenhonr excitation of the limit cycle when roll hysteresis 
Was applied. Yaw hysteresis, on the other hand, acted directly 
through the directional stability derivative, Vg, to produce 
a yawing moment. 

Seounmary of Hysteresis Induced Wing Rock 
There are several significant points worth considering 
from this analysis and the data presented in Table I and 
Picanese 59. These anclude: 


(1) The amplitudes of the limit cycle were directly dependent 
upon the magnitude of applied hysteresis. 


-(2) The period of the limit cycle was independent of the 
magnitude of applied hysteresis. 
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(3) While the roll rate was 180 degrees out of phase with 
pme sideslip angle tor the Dutch-Roll oscillation, it 
Piimeceerly INepiascmsor the limit cycle due to roll 
bys tenes s . 
fol acumeron to Piase Gditterences, the ratio of bank 
angle to sideslip angle was three times greater for 
the limit cycle due to roll hysteresis than for the 
Deven Roll oseillacion. 
moo lhewdaimit cy¢le for yaw hysteresis is essentially a 
forced Dutch-Roll response. Frequency and phase angle 
relationships are unchanged from the Dutch-Roll. 
eyimeneversing the relay action dependence on the sign of 
B results in a deadbeat response. Therefore, a wing 
rock limit cycle is dependent upon the sign of roll 
Or yaw hysteresis. 
(7) Combined yaw and roll hysteresis produces a limit 
cycle response which is essentially a compromise or 
average of the two effects. 
C. ANOTHER FORM OF AERODYNAMIC NONLINEARITY 
1. Previous Research 

As stated earlier, some research has been done and 
reported by Ross [Ref. 3] on wing rock as a limit cycle caused 
by nonlinear aerodynamics interacting with an unstable Dutch- 
Roll oscillation. Ross reports that a wing rock limit cycle 
occurred on the Handley Page 115 testbed aircraft as the 
result of increasing the angle of attack to a critical point 
where any small lateral-directional perturbation excited an 
alternating wing rock motion, which steadily grew to a con- 
Stant amplitude with the controls fixed. The motion would 
immediately damp to zero when the pilot commanded a reduction 
in the angle of attack below critical. The Handley Page 115 


is a British aircraft designed to serve as a flying testbed 


primarily for research on the Supersonic Transport and has 
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Meemeeric Characteristics similar to many current tactical 
@ircratt, i.e,, thin swept delta wings and a slender fuselage. 

In modeling the motion observed in flight and wind- 
mummel tests, Ross essentially took the linearized equations 
describing lateral-directional motion and imposed a nonlin- 
earity on yawing moment due to sideslip, Nw. The extent of 
eie nonlinearity was reflected in a polynomial of degree 
SUfeacient to match flight and wind-tunnel data which indicated 
a nonlinear variation of aerodynamic moments and forces as a 
memeeton Of Sideslip angle. A fifth order polynomial (in odd 
powers only) was proposed initially; however, as Ross pointed 
Siemmea cCUDIC nonlinearity was sufficient to yield good agree- 
ment with experimental results of wing rock. Ross concedes 
that although nonlinearities existed in the dihedral effect 
and sideforce characteristic, these were considered relatively 
insignificant and thus ignored. 

An analysis follows which resulted from an intention 
to simplify the mathematical modeling proposed by Ross and 
provide a better analytical tool for use in analyzing lateral- 
directional motion at angles of attack near stall. As the 
reader will see, the analysis and conclusions presented here 
do not necessarily agree with those of Ref. [3]. 

2. Theory 

As a means of simplification, the sideforce equation 
mcmenle accompanying Spiral divergence mode were, as in the 
hysteresis evaluation, assumed to be insignificant and thus 
eliminated from the equation set describing lateral-directional 
motion. This is not unreasonable, as flight tests involving 


45 





Sererart exhibiting wing rock have shown that the sideforce 
Paeietveteney Ssialleas to be non-discernible to the pilot 
Maaeheas virtually no effect on the heading or ground track. 
Therefore, the equations of motion applicable to this analysis 


Meeetic same as in the section on hysteresis, 1.e€., equations 


be 


b= -N, 8 ts -w%P ( 03a) 
y 


Ross concluded that the term which contributed most 


J 
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Lae — Lj G op ( 03b) 


to the limit cycle in view of an unstable Dutch-Roll oscilla- 
tion was a cubic nonlinearity in yawing moment due to sideslip. 


Pommenis analysis the cubic nonlinearity was introduced as: 
Wewing moment due to sideslip = Ng & = Ng (/# EB) 


where € is a scalar of arbitrary magnitude which facilitated 
varying the desired level of nonlinearity for a given analysis. 
[It was necessary to reduce the set of simultaneous differential 
equations in (05) to a single third order equation as a func- 
tion of a single independent Vecieaic Ole rao0 From equation (03a) 
Bicmclito le nonlinearity yields: 

BS =-NeB(MKEB) EN, 8 >No P (03c) 
ieee sence Gecrivative Of this equation with respect to time 


yields: 


(12) 
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Solving for / in (12) aswarpreeursOa tO Climination of the 


variable # SEOMmuncmeGuatuwens set yields: 


p 5 ig [-8 “Mg (/ 4 FEBS +N, ef (13) 
Now solve for # from equation (03a) and substitute the result 
in equation (03b). 
fhe bh br fe thgP 
2 4gBh bot LB -ng (esata) 04 


Equations (13) and (14) can now be equated to complete the 
rrnacton Of gvand yield a single third order differential 
equation in the single variable f/, from which the analysis 


proceeds. 
wey ‘ a : ys 4 . 
lB Ng (1436 BY) 8 +N, B/hgO-L, B — [:8-Ng (68 )8+N, 4] 
@learing through by a minus Vy Piece ecctine Cerms yreids. 
B-(N, +L) +[Ng (1+ 2€ 7°) “Nok, tN, Lol é + Mehg-4pNg Cr+ (B28 +0 


It is convenient to introduce the following notation for ease 


Om writing: 


A = = Ne + LJ where: b Ng vee es INT 5 


1 Yr-p pr 
B= b, f+ bp b, = €Ng 
ee Cap Ey = Nl aoe - LN g 

Cz = ~ENgL,, 


Thus, the resulting third order equation can be expressed in 


peuone general form as; 


3 Z 
of Tf + HE + oa) +2 ont 





This is a nonlinear differential equation for which there is 
MO analytical solution. Necessary assumptions were made in 
order to gain an approximate solution as discussed in the 


following. 


- 


Ss. Analysis 


With the resulting general form of the equation of 
motion given by (15) it was possible to proceed with an analysis 
of the motion similar to that presented by Ross, differing only 
fMmene degree of simplicity. Recall that the simplicity was 
gained by making the assumption that ¥= -Mor that the air- 
craft remains on a steady track (without heading deviation) 
while engaged in wing rock. Ross, however, did not make this 
assumption but rather included the side force equation result- 
ing in a fourth order system. It was further assumed, as did 
Ross, that the oscillatory motion regardless of damping could 


Demeepresemted by the following: 


i) 


iG ot a wt =VT cos 8 
peck — Vey = m~”/ 7 oR gs AT (16) 


Oya Ww 


| 


where v 


i) 


© r=) 


This simply says that wing rock motion can be described as a 
Sinusoidal oscillation with sideslip amplitude considered as 
a Slowly varying function of damping and time. The shorthand 
notation of sigma (c) and theta(@) was used for writing 
SGmvenience. 

In order to gain an approximate solution to equation 


(15), the Krylof and Bogoliubof (KB) method, Ref. [7] was used. 


48 


67 


7 © 
ci < <4)! 1” 


° Vs ta! 


114 


& 2.64 





iias method Nas been used extensively to analyze second order 
oeeems Witn “weak” nonlinearities. With the exception of 
Mies analysis, as well as that done by Ross on the fourth 
Meet eoyocetr, NO Other reference to use on higher than second 
order systems could be found. Nevertheless, there was no 
apparent mathematical restriction to applying the KB method 
to the third order system of equation (15). Moreover, the 
KB method is applicable given the assumption that the ampli- 
tudes of a vibratory system vary only slightly over a single 
Sete. Inus, the terminology often used is the assumption of 
Slowly varying parameters, which is equally valid when con- 
Sidering a slowly divergent Dutch-Roll oscillation. 

It was necessary to take the assumed form of the 
solution in equation (16) for the dependent variable, @, and 


derive the first, second and third derivatives, which were: 





= 97 [A cos® SR sag) | (17a) 
2 
va =¢ [(Atut) cose - 24Asine] (17b) 


rT [ (A=3W*A) cosd + (w= FWA") sine] (170) 


the 
eh 
TT 


PMibcelemting equations (17) into equation (15) yields: 


: (*f dB) 
— 2 Sf. EL. ce) 





Sat asdfS [reo oen 0 


2 we 
-34,] A C05 @- WSN e/r cbs 8~ ¢, Case-C;,T Cos of 


de? rife web Lee, cos + [eawd + b,c J sm s 


Q. 

ft 
w 

il 


-134, 4 #0, |r *cos%e + 3b, Lat sx & Cos "@ f (18) 


49 





memetion (18) is of the form: 
3 
of = TFA w, T, cos 9, syv 6) 

Mme eoverning principle of the KB method says that since the 
amplitudes are considered essentially constant over a given 
mele (slowly varying parameters), then the average value 
over the cycle is a sufficiently close approximation to the 
actual maximum value. This averaging process, when applied 


meomeguation (18), yields the following: 
S z 
[7-30] 
[suk] = #["" 
Ww" Puy A 7) Fsvbde 


[U?- sul] = Bb-T[ak* ut) th dee] Zrrfstedtq]e* 
(19a) 
fw i sui] 2 o) [eau x b wo ) + trbwr'} (19b) 


Arriving at equations (19) necessarily involved the orthogon- 


277 
f 
z/ F tos ade 


il 


aimeey principles of integration as shown below: 
27 27 
f Cos 6 de I siy* @ oa = 77 


27 227 
4 2 
yA. do > ta] siv 8 te88 d&8 50 
0 ro) 


2/7 27 i 

3 t TP 

bi Siv8 (os°@ A820 ANG aa cas "6 do = yn) 
Note that equations (19) yield the characteristic equation 
for the stable system, when the nonlinearities vanish which 


abe contained in the odd subscripted variables, bz and Cz: 
Sy 3 ‘ 
(Aviw) = LA? 3u%hf- iw ss? ] 
S; , 
aaa (x tlw) 2-4 (Lt =, (d-eiu)- ¢, 
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which is equivalent to: 


3 
So eae ee en 


Peoumameethat a limit cycle occurs, the following conditions 


must hold: 
yee Te SE te, 
ea 
G = De 
Thus Ine = 0 = damping 
7; = D = constant 


where the subscript, s, is used to denote steady-state condi- 
tions. Hence, from equations (19), with the limit cycle 
conditions of zero damping imposed, As = 0, there is 

ais 


Oi Als ~ €, - % C3 Ts 


a 


F. G 
Ww = bw +t “&S, WO; 


ftemeda be reasonably assumed that a limit cycle will not occur 
for a frequency “= 0. The phrase "reasonably assumed" is used 
gameesedly. Although it seems intuitive that such an assumption 
would be valid, some researchers, notably Parkinson [Ref. 8] 
detailed a study of the transverse galloping of a long prism 

of square section in a steady flow field in which the subject 
body possessed an unstable focus at rest. Whether this means 
ioe os a Valid solution for a wing rock limit cycle, 
Memeitinecannot be discarded, remains largely unanswered. 
Nevertheless, such an assumption was made, as did Ross, enabling 
elimination of the & term between the two above equations 


Weta? tle Lindl Caquations describing the limit cycle. 


Spb 





2 J = 
Alb, + bet )- G&G = 





z 
Co, GP, + Ae, - ©, = © 
- Ab, -%, 
2 I le: = 
/ 
Ab, - 5 a 
v5 = pice (Y=) (20) 


Meemerequency, ©, tor the limit cycle can then be obtained 


as. / 
4b -S ie 
aye (he tbs (EE) a 


Meee that one of the preconditions for the limit cycle as 
proposed by Ross was that the aircraft must possess an unstable 
Dutch-Roll mode of slowly varying parameters, which is con- 
strained in amplitude by a positive cubic nonlinearity in 
yawing moment (yawing moment increasing with sideslip angle). 
Vemmenere was nothing explicitly incorporated into the mathe- 
matical model, which yielded equations (20) and (21), that 
would account for this condition. Instead, it was assumed 
thatea limit cyche existed due to a constrained motion in- 
Stability. The modeling proceeded on this basis with the 
Metron Comstraint explicitly incorporated by Ng (tr 6 )6 . 
Tieemonoula not pe interpreted to mean that a limit cycle 
will occur only if there is a divergence of the vibrating 
System initially (refer to section III.A); however, when 
applied to aircraft dynamics this is probably the most common 
Tyee cyene Oscillations occur, as was experienced by 


nh@ssmancde nas been experienced by the author as well. 


Sf 





The effect of inducing a limit cycle via a positive 
nonlinearity on Ng at high angle of attack is intuitively 
aepealing. This would amount to increasing directional sta- 
Dableet ye with Sideslip, or the spring in the spring mass damper 
Meee stitfer with displacement from equilibrium, until there 
is an equilibrium established between the energy contained in 
tie Unstable Dutch-Roll and that absorbed by the spring. 
Pimestedlly, this could correspond to an aircraft exhibiting 
a decrease in directional stability at low sideslip angles 
and high angle of attack due to a turbulent, low energy flow 
field emanating from the wing and engulfing ehe empennage 
section, thus reducing the restoring moments generated by 
Seo@olrzing and/or control surfaces. This is often the case 
fOr the tactical high performance aircraft. However, as the 
tail is yawed out of the low energy flow and into a lower- 
timembence-level, higher energy flow (higher kinetic energy 
or dynamic pressure), the directional stability will increase. 
Thus, the premise of an unstable motion constrained by a non- 
mitedmity in directional stability would appear to be valid. 

The F-94 airplane as originally used had a stable 
Dutch-Roll; therefore, it was necessary to destabilize the 
airplane artificially for the sake of conducting wing rock 
studies. The primary means of achieving an unstable Dutch- 
Roll oscillation consisted of altering the sign and magnitude 
of the yaw rate damping derivative, N... The effect of this 
Can be seen in figure 10, where the motion slowly diverges at 


a frequency equal to the stable Dutch-Roll. Since the 
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denominator of equation (20) is positive, and in order to 
ensure that the radical term remained positive, it was neces- 
fey to alter a second derivative such that Ab, > Cy: After 
some manipulation, the following derivatives were changed as 
mmicdacated: 

Nee. OZ 


Np a te 


a. Predictions Based on Analytical Solution 
The following 1S a sample calculation of a pre- 


Gemeced limit cycle using equations (20) and (21). 


> 


= -( Neer L =e eS) 
(Ny + Lp) 


Nee Nei = lene 


Q 
I 


i ee = Nel, = 4.0916 


oO 
li 


ENg Sly Sila 


3 
Cz =  €Nel,, Se els y, 
ry = ase (Gass OsT ety 0948) 
with € = 5 

Me = £1267 Rad 47.15 Deg 

24 =e 5454 Rad/sec 

is = 490 / SCC. 


In addition to the analog and digital computer 
programs previously mentioned, programs were likewise written 


and executed based upon the reduced third order equation (15). 
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RADIANS 


These are presented in Appendix D and F. These programs were 
mMsed primarily to compare with the results based on the two 


-aeotton System (equations 03b, Q3c) to ensure consistency 


ame accuracy to the maximum extent possible of the mathematical 


megeling that led to equations (20) and (21). The same plot- 
ting program of Appendix C was used. Heavy reliance was 
placed upon the analog computer for qualitative analysis, 
where changes could be made and their effects observed imme- 
Gmately. 
Doe CMemecatTon or Predictions 

The limit cycle predicted by the above sample 
Calculation was not obtained. Rather, another limit cycle 
of lower steady state sideslip amplitude, figure 10, occurred 
almost immediately upon release from initial conditions with 


a level of nonlinearity, €=5. 


UNS TALE (DIVERGENT) DuTeH=@akd 
—-~—~ CONVERGING LIMIT CYCLE 


a e BS / Fe ay om i 1/2 
/ { J \ ey 
J P / f \ } \ 
J ‘ f f f i , 
2.gn ! f 
\ ; \ \ 
\ \ es 
-2.25 \ \ \ t \ f 1 / 
a Mo WW te / 


rm is cS x = 
TIME (SEC) 


ooo Wire History of Wing Rock Limit Cycle 
and Divergent Dutch-Roll. 
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The value of the predicted go = .13 radians was greater than 
the initial condition fy = -08. Thus, according to theory and 
the fact that an unstable Dutch-Roll existed as shown in Fig. 
10, the motion should have slowly diverged to the predicted 
Steady state values. Similarly, an initial condition greater 
team the predicted value of Te was tried (not shown) to see 
af the motion would converge to the limit cycle; however, the 
time history was essentially identical to the first in that a 
limit cycle again developed with an amplitude very nearly the 
Same as the initial conditions. Therefore, the wing rock 
Moe cycle appeared to be a function of the stability deriv- 
and N anda taeceintedal €ondrtions, rather than of 


is (eh 


the value of the nonlinearity in yawing moment. 


atives, N 


The effect of the level of nonlinearity as viewed 
on the analog computer was merely to alter the frequency, with 
no change in amplitude. On the other hand, the predicted 
frequency, We» agreed very well with experimental data. For 
the above sample, the predicted WwW. in terms of the period 
(where period We = a was 4.67 sec, While experimental results 
Showed 4.7 sec. This Close agreement on frequency held for 
changes in the level of nonlinearity €&, as well as changes in 
the stability derivatives, as long as the requirement Ab, > Cy 
was satisfied, as previously discussed. 

Experimentation on the analog computer showed 
Mite omemc emus tony twesponse behaved as if neutral stability 


existed. When considering motion without nonlinearity imposed, 


the Routh Stability Criterion can be used to determine the 
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Hemera stability boundary, and it appeared at first that 
Peaetbly ene Rowth discriminant had been inadvertently satis- 
fied through the manipulation of Ny; and Ne? me. , Ab, - Cc, = 0, 
i directly proportional to Routh discriminant. However, this 


was not the case as Ab, - cy was) foundmequal tomo. 1158. 


4. Summary of Wing Rock Induced by Nonlinear 


Yawing Moment 

Now, an obvious question arises as to the source of 
disagreement between these findings and those of Ross. The 
main difference with this analysis and that carried out by 
Ross was the simplification imposed at the beginning, i.e., 
that the sideforee equation and accompanying spiral mode were 
insignificant and need not be considered. Since the F-94 
possessed no natural instability in the Dutch-Roll at the 
angle of attack tested, the artificial manipulation of the 
Seamrlity derivatives N, and No? as discussed above, may have 
been sO arbitrary as to not correctly represent the character- 
istics previously exhibited by the Handley Page 115, and docu- 
mented by Ross. Whether these or other concepts were the 
source for the disagreement in results is unknown. Possibly 
the analysis Ross presented for the Handley Page 115 was 
peculiar to that aircraft, although intuitively one would 
tend to argue in favor of the premise that a limit cycle 
would occur for an unstable Dutch-Roll in the presence of an 
increasing yawing moment with sideslip angle. In view of 
this, it is recommended that follow-on research be conducted 


Using the complete set of equations of motion, i.e., including 


oy 





the sideforce equation and the analysis technique presented 
here to make a more conclusive determination as to the validity 
O@ the concept that a wing rock limit cycle can be the result 
Seen unstable DutchwRoll coupled with a cubic nonlinearity 

in Ng. Moreover, any additional research should be conducted 
On an aircraft known to exhibit wing rock motion such as the 
hemabey Page 115 or preferably on an operational, tactical 


@mperatt such as the Northrop 1-38 or F-5. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


The existence of a lateral-directional aircraft limit- 
meee oscillation, due to an aerodynamic hysteresis corres- 
Pomeeme to a roll or yaw moment rebay action that is dependent 
eon the sign of sideslip velocity, has been demonstrated. 

The use of state variable notation allowed the limit-cycle 
motion to be expressed in a concise and clear manner. 

Although bank angle was not a state variable in the 
@melysis and hence could only be estimated from roll rate, 

@he inclusion of this term by the introduction of a sideforce 
equilibrium equation is practicable for a more detailed anal- 
ysis. However, in keeping with the purpose of the analyses 
presented here, to illustrate a concept, the only equations 
used were those absolutely required to achieve that goal. 

The existence of a wing rock limit cycle oscillation due 
to a nonlinear yawing moment, as reported by Ross [Ref. 3], 
was not confirmed by the analysis presented here. The essence 
of the differences in the two analyses was in the basic assump- 
tion that the sideforce equilibrium equation need not be 
included. Such an assumption is based in large part upon 
Pena@nalerlieht test experience of the author, in which the 
Characteristics of wing rock at high angle of attack have 
Clieagkyernolcated that yaw excursions are relatively small, 
thus adding support to the contention that sideforces be 


assumed negligible. 
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From the exhaustive experimental observations made in 
the laboratory, albeit using an aircraft that was known not 
memDOSSesSs any Wing rock limit cycle characteristics, serious 
doubts have been raised concerning the generality of Ross' 
omeeiysis. While on the one hand the concept of a nonlinear 
poems, aS applied to aircraft dynamics, and the "limiting" 
@eechne divergent Dutch-Roll oscillation, has intuitive appeal, 
Experimental results using analog or digital modeling consist- 
Smey Indicated an effect only on the frequency of the vibra- 
Beene With no effect on damping. Indeed, a limit cycle was 
achieved based upon the modeling criteria used by Ross, but 
iemmase predicted with respect to steady state amplitude. 
However, the predicted and resulting experimental value for 
frequency of the oscillation showed excellent agreement. Of 
the limit cycle obtained, the damping appeared to be a func- 
tion of the modified stability derivatives necessary to 
artificially destabilize the Dutch-Roll of an otherwise 
stable test aircraft (F-94) rather than of the level of the 
Subic nonlinearity. 

iGessOusce Of diSaemecment between the results obtained 
Meaemana those Of ROSS remained unknown at the termination 
of this research. Sufficient time was unavailable to pursue 
the question further; however, it is recommended that follow- 
on research be conducted, when time permits, to make a more 
deedileawandlysis, possibly including the sideforce equa- 
BIoieetimanecrrOrt to determine the source of disagreement 
between the two analyses, which has led to quite different 
aircraft responses. 
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we Lewis enOpm@deciac these analyses will make clearer 
tMat interpreting an actual aircraft time history as represent- 
ing a Dutch-Roll mode with exactly zero damping may be mis- 
Peiding, Since a mOtion with very similar appearance could be 
memmmt Cycle due to aerodynamic relay, or hysteresis actions, 
Or nonlinear yawing and rolling moments as posed by Ross. The 
mimeeect problem of determining the time response of an aircraft 
iemeo LOrcing functions, be they linear or nonlinear, is 
Meeecaple. The inverse problem of identifying the plant from 
the aircraft time histories (known as aircraft parameter iden- 
tification) is much more difficult, because of the uncertainties 


Mamestilidting types of plant nonlinearities. 
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APPENDIX E 


ANALOG COMPUTER PROGRAM FOR 


-Ng (1+ cB°)B +N,6 - NP + AN*SIGN (8) 


es 


bob adh ise Lyp + AL * SIGN(8) 


p 


gs7ogwis $9 NOLLINSYIO HOS + XIONFAAV SIS 


J7 





70 





ae eee X. F 


ANALOG COMPUTER PROGRAM FOR 


SL1TOA O] SNNIW YO S/d (+) 


YOLYYDILNT <|} 


(WS ANI) 89171 78WY ra: 
YILIWAOILN I LOS C) 


ye eS "TINW GILYIANI | 1 | 
? YITTIATLINWG [| 


Q/ 





B= -~AB - 8 (by + 3b38°) = ¢1B- c3B> 


dyénC ~ t) SRWIS 
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